Non-relativistic limit in the 2+1 Dirac Oscillator: A Ramsey Interferometry Effect 
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We study the non-relativistic limit of a paradigmatic model in Relativistic Quantum Mechanics, 
the two-dimensional Dirac oscillator. Remarkably, we find a novel kind of Zitterbewegung which 
persists in this non-relativistic regime, and leads to an observable deformation of the particle or- 
bit. This effect can be interpreted in terms of a Ramsey Interferometric phenomenon, allowing 
an insightful connection between Relativistic Quantum Mechanics and Quantum Optics. Further- 
more, subsequent corrections to the non-relativistic limit, which account for the usual spin-orbit 
Zitterbewegung, can be neatly understood in terms of a Mach-Zehnder interferometer. 

PACS numbers: 42.50.Vk, 42.50.Pq, 03.65.Pm 



INTRODUCTION 



The natural relativistic extension of the quantum har- 
monic oscillator, known as the Dirac oscillator has 
become a cornerstone in Relativistic Quantum Mechan- 
ics. It was initially introduced as a relativistic effective 
model to describe mesons, since it presents interesting 
quark-confinement properties 0, Q. Moreover, subse- 
quent studies have revealed several amazing properties of 
the Dirac oscillator in different contexts. Beyond its ex- 
act solvability, the energy spectrum presents certain pe- 
culiar degeneracies which can be related to a non-trivial 
symmetry Lie algebra Furthermore, its solvability 
can be traced back to an exact Foldy-Wouthuysen trans- 
formation 0, and its special properties are related to 
a hidden supersymmetry Q. Additionally, the positive- 
and negative-energy solutions are associated to super- 
symmetrical partners, which ensures the stability of the 
Dirac sea under the Dirac oscillator coupling 0]. 

Some analogies between the dynamics of this relativis- 
tic model and the typical Jaynes-Cummings (JC) dynam- 
ics in Quantum Optics @ have been discussed in [3| , and 
some of its non-relativistic properties in 
markably, in a two-dimensional setup this analogy be- 
comes an exact equivalence between the 2+1 Dirac oscil- 
lator Hamiltonian and the Anti-Jaynes-Cummings (AJC) 
interaction [l^ , which builds a bridge between two unre- 
lated fields. Quantum Optics and Relativistic Quantum 
Mechanics, and favors a fruitful exchange of ideas be- 
tween both communities [3] . Relativistic effects such as 
the Zitterbewegung, a helicoidal motion performed by the 
average position of a free relativistic fermion, can be rein- 
terpreted with the language of Quantum Optics. This 
dynamical phenomenon also becomes observable in the 
spin and orbital degrees of freedom, where it can be sur- 
prisingly interpreted in terms of optical Rabi oscillations. 
In this work we shall be concerned with the extension 
of this novel perspective onto the non-relativistic limit 
of the two-dimensional Dirac oscillator. In this manner, 
we are able to identify a novel feature of the Zitterbewe- 
Qunq, which is interpreted as a Ramsey interferometry ef- 
fect lia . This remarkable effect shows that the standard 



non-relativistic limit described in Relativistic Quantum 
Mechanics textbooks should be reconsidered in this un- 
usual scenario p^ . 

This paper is organized as follows: in Sect. |TT1 we re- 
view the properties of the two-dimensional Dirac oscil- 
lator and the exact mapping onto an AJC Hamiltonian. 
In Sect, nil] the non-relativistic limit is considered from a 
Quantum Optics perspective, which allows the prediction 
of a novel kind of Zitterbewegung which is interpreted as 
a Ramsey interferometric phenomenon in Sect. IIVI In 
this section we also discuss how the effects of this inter- 
ference process have strong consequences in the electron 
trajectory. In Sect. |Vl the analysis of additional correc- 
tions to the non-relativistic limit is discussed in terms of 
a Mach-Zehnder interferometer. There, we find that the 
first order non-trivial correction already shows a pertur- 
bative spin-orbit Zitterbewegung Finally, we con- 

clude reviewing the consequences of our work in Sect I VII 
In appendices and |B] we give detailed derivations of 
the standard non-relativistic limit and its complete per- 
turbative series, respectively. 



II. TWO-DIMENSIONAL DIRAC OSCILLATOR 

The physical laws that describe the properties of mi- 
croscopic particles are accurately described by Quantum 
Mechanics, and in particular by the Schrodinger equa- 
tion. Nonetheless, quantum phenomena occurring at 
high energies cannot be properly addressed by such the- 
ory, and one must employ Relativistic Quantum Mechan- 
ics jl^. A relativistic spin-1/2 particle of mass m is de- 
scribed by the Dirac equation 



[ca • p + /3mc^] |*), 



(1) 



where \^) is the four-component Dirac spinor, Ui := 
off-diag(cri, cTi) , and (3 := diag(l2, — 12) arc known as the 
Dirac matrices which can be expressed in terms of the 
usual Pauli matrices ct^, p is the momentum operator, 
and c stands for the speed of light. 

The Dirac oscillator is obtained after the introduction 
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of a peculiar coupling in the above equation ([T|) 



dt 



(p - im/3wr) + /3mc^] |*), (2) 



where u stands for the Dirac oscillator frequency, and r 
represents the particle position. The relativistic coupling 
p — im/Swr, known Dirac string, cannot be under- 
stood as a simple minimal coupling procedure, and is re- 
sponsible for the special properties of this relativistic sys- 
tem. In particular, the non-relativistic limit of the afore- 
mentioned Dirac oscillator ^ leads to the usual non- 
relativistic harmonic oscillator with an additional spin- 
orbit coupling which shows the intrinsic spin structure of 
the relativistic theory 

The restriction to two-spatial dimensions appreciably 
simplifies the relativistic problem, since the Dirac matri- 
ces become 2x2 matrices which can be identified with 
the so-called Pauli matrices = cFxiCty = Cy, (3 = a^. 
In this manner, |^) can be described by a 2-component 
spinor, and the Dirac oscillator model now takes the form 



dt 






l^>- (3) 



This two-dimensional system was algebraically solved 
by introducing chiral creation and annihilation op- 
erators 



ai 



where ax,al.,ay,a^y, are the usual annihilation-creation 

i-^p-'), and A = ^Jhjmljj rep- 



operators a 



71 Va' 



resents the ground state oscillator width. These opera- 
tors allow an insightful derivation of the energy spectrum 



(5) 



where the integer ni stands for the number of left-handed 
orbital quanta, and ^ := hw/mc^ is an important param- 
eter that specifies the importance of relativistic effects in 
the Dirac oscillator. A different approach, based on the 
solution of differential equations, has also been discussed 
in .The associated eigenstates are found to be 



(6) 



where \x^) ■= (1,0)* and |xi) 
the Pauli spinors, and a„, {Em + mc'^)/2Eni and 

St ants. 



(0, 1)* are known as 



mc^)/2Eni are real normalization con- 



The notation in Eq. ^ clearly shows that the rel- 
ativistic eigenstates exhibit entanglement between the 
orbital and spin degrees of freedom. This entangle- 
ment property is essential in order to obtain the rel- 
ativistic effect of spin-orbit Zitterbewegung, where cer- 
tain oscillations in the orbital and spin angular momen- 
tum arc unambiguously identified with the interference 



of positive- and negative-energy components. Introduc- 
ing the spin and angular momentum operators Sz = ^ha^ 

and Lz ~ h{alar — ajai), and considering the initial state 
l'I'(O)) jn; — l)|x|), one immediately finds the so-called 
spin-orbit oscillations associated to the Zitterbewegung 



iL-)t - -TTl|77^sin2c^„,t - h{ni - 1), 

{Jz)t^ 



(7) 



where Jz ~ Lz + Sz stands for the z-component of the 
total angular momentum and is obviously a conserved 
quantity. 

Finally, we recall the interesting mapping between the 
two-dimensional Dirac oscillator and the A JC model [3] , 
where the relativistic hamiltonian can be written as 



H = h{g(7^aj + g*(T a;) + 



6a, 



(8) 



with cr+ := |xT)(Xi I, cr -.^ \xi){x]\ ss the spin raising 
and lowering operators, g := 2imc^y^/h as the coupling 
strength between orbital and spin degrees of freedom, 
and where 6 := mc? can be interpreted as a detuning pa- 
rameter. In Quantum Optics, this Hamiltonian describes 
an Anti-Jaynes-Cummings interaction, and can be im- 
plemented with trapped ions . Within this novel per- 
spective, the electron spin can be associated with a two- 
level atom, and the orbital circular quanta with the ion 
quanta of vibration, i.e., phonons. Note that there is also 
the possibility to map this relativistic Hamiltonian onto 
the more standard Jaynes-Cummings model [3. In the 
following section, wc derive the non-rclativistic limit of 
the two-dimensional Dirac oscillator ([5]), and discuss the 
nature of the physical properties described in Eqs. ©-([T]) 
in this non-relativistic scenario. 



III. NON-RELATIVISTIC LIMIT IN QUANTUM 
OPTICS 

The original Quantum Optics perspective of the two- 
dimensional Dirac oscillator in Eq. ^ stimulates the use 
of quantum optical tools in a relativistic quantum frame- 
work, and viceversa. In particular, we can use the quasi- 
degenerate perturbation theory [l9| in order to derive an 
effective Hamiltonian in the non-relativistic limit. This 
regime is attained when the relativistic parameter fulfills 
^n; <C 1, which allows the usual description of the Dirac 
oscillator Hamiltonian in Eq. ([H]) 

H = Ho + XV, (9) 

where Hq is the unperturbed Hamiltonian 

Ho = Saz, (10) 

and XV represents the following perturbation 

XV ^ X{a+al -a-ai), (11) 
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where the interaction coupling A := 2imc^v^ satisfies 
|A| ^ (5, and consequently the coupling in Eq. ([TT|) can 
be treated as a small perturbation. In this regime, the 
Hilbert space can be divided into two approximately dis- 
connected subspaces 



which can be readily evaluated as 



where 



W| = span{|r^^)|xT);?^; = 0,1...}, 
Hi = span{|n;)|xj.);n; = 0,1...}. 



(12) 



(13) 



This is easily understood from the fact that ^ni <^ 1 
is equivalent to huj S, which implies that the pertur- 
bation in Eq. pT|) does not suffice to induce transitions 
between the spinor components. In this manner, the sub- 
spaces corresponding to the spinor degrees of freedom 
become decoupled, and we may describe the effective dy- 
namics in such subspaces ( see fig.[I]). 





lAI <(5 



FIG. 1: "(Color online)" Schematic diagram of the Dirac os- 
cillator couplings in the non-relativistic limit. The subspaces 
Til ^'id 7Y| become disconnected in this regime 



2hujajai 




(16) 



The effect of the relativistic corrections in this regime 
can be understood as a level shift with respect to the rest 
mass energy that depends on the number of left-handed 
quanta. The non-relativistic energies associated to the 
corresponding eigenstates |»t.z)|xt) ^^^^ I"' ~ l)IXi) 



^I"I>|XT> = +'71C^(1 + 2^"i)i 



(17) 



which are equivalent to the leading order correction of 
the exact eigenvalues ^ in the limit ^ni <C 1, namely. 



E = ±mcVl + 4^^/ ~ ±mc^{l + 2^ni) + O ((C";)^) • 

(18) 

Therefore, we obtain the so-called energy shift term 
AE = 2mc^(l -|- 2^71;) , which is usually known as a 
dynamical Stark shift term in the quantum optics litera- 
ture [1^. In Optics, the non-relativistic effective Hamil- 
tonian (jl6p achieved for large enough detuning is equiv- 
alent to the dispersive linear susceptibility and the real 
part of the refraction index, with opposite contributions 
from the excited and groimd states [21| . 

This noteworthy interpretation of the non-relativistic 
limit in terms of measurable optical quantities is shown 
to be equivalent to the standard non-relativistic limit in 
Relativistic Quantum Mechanics in Appendix [XI In the 
following section, we shall use this remarkable perspec- 
tive to describe a novel sign of the Zitterbewegung, which 
can be understood in terms of a Ramsey interferometry 
effect fii. 



IV. ZITTERBEWEGUNG IN THE 
NON-RELATIVISTIC LIMIT 



In order to obtain the effective Hamiltonian, we rewrite 
the perturbation in Eq. ([TT|) as follows 



(14) 



where fi = 1,2, the operators Ai = ai, A2 = a\ describe 
the slow varying orbital degrees of freedom, whereas the 
operators Bi = a~ , B2 = cr"*" represent the coupling of 
the fast spinorial degrees of freedom, and Ai = g, A2 = 5*. 
The effective Hamiltonians correspond to 

His^-S + J2{xi\X^B^j^^\^'B^,\xi)A^A^,, 



(15) 



As discussed in previous sections, the interference be- 
tween positive- and negative-energy components gives 
rise to a relativistic oscillatory behavior known as Zit- 
terbewegung. This phenomenon has a pure relativistic 
nature, and therefore it is usually believed to vanish in 
the non-relativistic limit. Contrary to common belief, 
we show in this section how a peculiar Zitterbewegung 
effect can still arise in the non-relativistic regime if the 
initial state is appropriately prepared. Furthermore, we 
also discuss how this dynamics might be interpreted as 
a Ramsey interferometry phenomenon, and how it can 
lead to measurable effects in the particle trajectory. 

Let us consider an initial state |5'(0)) := Q^|?^^)|XT) + 
I3\ni — l)|xi): which involves both spinorial components, 
where a, /? G C are correctly normalized \a\^ + |/3p = 1. 
One directly observes that this state mixes the positive- 
and negative-energy solutions in Eq. ()17p . which is the 
fundamental ingredient leading to the Zitterbewegung. In 
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order to obtain such effect, we derive the time evolution 
under the effective Hamiltonian 

|vl/(t)) =«e"'^^",*|„,)|;^^)+/3e+'^",*|n,-l)|xx), (19) 

where 57„, :— rac^(l + 2^ni)/h. The corresponding spin- 
orbit expected values ([7]) become 

(L,)t = -ft(n, + |/3|2), 

{S.)t = I (lap - m , (20) 
(J.)* = f (2n, + |a|2-3|/3|2), 

where any remainder of the original oscillatory Zitterbe- 
wegung in Eq. ([7]) has completely vanished. Nevertheless, 
positive- and negative-energy components are simultane- 
ously involved in the initial state, and therefore there 
must exist some kind of novel Zitterbewegung. 

The possibility to observe such effect can be achieved 
if we consider the superposition of states with different 
energy modulus I-Enj- This novel Zitterbewegung can 
be exemplified by a Ramsey interferometer in which the 
field is prepared in a coherent state and the atom in a 
50 % superposition of its excited and ground states ( 
i.e. ajxt) + /3| y | )) by resonant interaction with a classi- 
cal laser beam [2]. Afterwards, the effective interaction 
produces a different field evolution conditioned to the 
atomic level. Finally, resonant interaction of the atom 
with a second laser beam mixes the contributions from 
the upper and lower levels leading to the interference of 
the positive- and negative-energy evolutions. Therefore, 
the Ramsey fringes such as the ones in Refs. [22| can be 
regarded as suitable examples of Zitterbewegung in this 
particular regime. 

Let us exemplify this Hamscy- Zitterbewegung with 
the following initial state prepared in a superposition 
of the two-spinor states |5'(0)) := \zl){a\x^) + /3|xj.)), 
where \zi) := e"'^'' X^n, / represents a left- 

handed coherent state with zi G C. This initial state also 
involves positive- and negative-energy components, and 
its time evolution is 

|*(i)) = ae-'"°*|z,e-2-*)|xT> +/3e+'^'^*|z;e+2-*)|xi). 

(21) 

As time elapses, the phase evolution of the orbital co- 
herent state is strongly correlated to the internal spino- 
rial degree of freedom, just as the Ramsey interferomet- 
ric time evolution. This peculiar correlated dynamics is 
a clear consequence of the coexistence of positive- and 
negative-energy modes in the initial state, and therefore 
it stands as a direct symptom of Zitterbewegung. The 
final step in the Ramsey interferometric experiment is 
to recombine both spinor components leading to the in- 
terference of positive- and negative-energy modes. The 
measurement of Sz after this mixing effect is equivalent 
to the measurement of the x-component of the spin an- 
gular momentum Sx '■= ^(Jx in the state (|2ip 

{Sx)t = ftV(t) cos [{no + Q,)t + \zi\^ smiAU) + arg(a*/?)] 

(22) 



where V{t) = |a*/3|e-2N'l' '^'"(2'^*) can be identified with 
a periodic visibility factor which precedes the desired 
Ramsey interference term. Therefore, an oscillatory be- 
haviour in the x— spin component can be directly traced 
back to a Zitterbewegung in the non-relativistic limit. 
Note that this phenomenon is consequence of an ap- 
propriate preparation of the initial state, involving both 
energy modes, rather than the pure dynamical effect in 
Eqs. ©. 

This Ramsey interferometric dynamics pip can lead 
to measurable effects in the electron trajectory, since 
the interference of both energy modes undeniably causes 
a deformation of the particle orbit. The electron tra- 
jectory is described by the following position operators 
X = A(aj. + al + ai + a\)/2, y = iA{ar - a}. — a/ + aj)/2, 
whose expectation values evolve according to 

{x)t = A\zi\ [\a\^ cos(2u;t - 4>i) + \(3\^ cos(2w< + cjji)] , 

{y)t = A\zi\ [\l3f sin(2cjt -f (j>i) - \af sin(2wt - (j>i)] , 

(23) 

where we have used z; = |z/|e"^'. The particle trajectory 
described in Eq. has a remarkable periodic charac- 
ter, and must be compared to the Zitterbewegung-free 
trajectory of an initial state |4'(0)) := |zi)|xt), which is 
described by a circular orbit 

{x)t ^ +A\zi\cos{2ujt~ (t)i), 
{y)t^-A\zi\sm{2Lut~(j)i). 

Comparing both trajectories in Eqs. we realize 

that the Zitterbewegung phenomenon leads to a deforma- 
tion of the electron circular orbit (|24p ( see fig. ^ 

In light of the results presented in this section, we may 
conclude that Zitterbewegung phenomena may also arise 
in a non-relativistic regime as long as the initial state 
involves both energy modes, which is still a relativis- 
tic property. The initial state, which can be described 
by a coherent superposition of positive- and negative- 
energy solutions, cannot be described in the realm of 
non-relativistic quantum mechanics. Therefore, the per- 
sistence of Zitterbewegung in the non-relativistic regime 
can be traced back to the relativistic nature of the initial 
state, when carefully prepared. 

Furthermore, we have also described how this interfer- 
ence effect can be unexpectedly interpreted in terms of 
Ramsey fringes in the context of Quantum Optics. The 
effect of this relativistic interference can be observed in 
oscillations in the spin component, or even more dras- 
tically in the deformation of the electron circular orbit 
into an elliptic trajectory. 

V. CORRECTIONS TO THE 
NON-RELATIVISTIC LIMIT 

In section IIIII we have discussed the non-relativistic 
limit of the two-dimensional Dirac oscillator, where the 
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decoupling of the spinor components leads to a dynami- 
cal Stark shift in the energy levels. It was precisely this 
energy shift, which allowed the surprising description of 
the Zitterbewegung in terms of a Ramsey interferometry 
effect in Sect. lIVI In the present section, we consider how 
relativistic effects modify such picture as the parameter 
^ increases, and interpret the usual spin-orbit Zitterbewe- 
gung described in Eqs. ([7]) as a first order perturbation 
term. The picture developed in section IIIII is no longer 
valid, since it assumes the decoupling of the spinor sub- 
spaces (|13p. which forbids this peculiar dynamics. There- 
fore, we discuss a novel description which allows the in- 
terpretation of the spin-orbit Zitterbewegung in terms of 
a Mach-Zehnder interferometer. 

Let us consider the AJC-mapping of the two- 
dimensional Dirac oscillator ([5|) , which allows the descrip- 
tion of the Hilbert space as a series of invariant subspaces 
^ = ©^^=1 T^ni , where 

-span{|nz)|xT),|n,-l)|xi)}. (25) 
The relativistic Hamiltonian in these subspaces reads 

Hni = mc^ ((72 - 7]ni(Jy) , (26) 

where we have introduced a parameter 77„j = 2\/^ni di- 
rectly related to the small relativistic parameter ^n; <C 1. 
This interaction can be considered as a rotation of the cr^ 
term along the x-axis 

H,,, = 1 + rjlc~''-'''^aj'"'^^, (27) 



where the rotation angle satisfies tan(26'„j) := rin^ [2^ . 
The unitary time evolution operator can be expressed as 
follows 

U^^ = c-^^"< = c-'^-i^^e-'"!'"' W'^'e''^"^"- , (28) 
where 

</>«,(*) = ^\/l + <i- (29) 

This evolution has a clear interferometric interpretation 
in terms of a Mach-Zehnder interferometer ( see fig. [3]). 
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FIG. 2: "(Color online)" Electron trajectory for different ini- 
tial states where the coherent state phase is (j>i — n/2, and 
ai = {1, with i = 1,2,3. Note that ai = 1 

corresponds to the circular Zitterbewegung-bee evolution in 
Eq. (|24|) . as compared to the elliptic orbits caused by a Zit- 
terbewegung interference. 



FIG. 3: "(Color online)" Mach-Zehnder interferometer dia- 
gram of the Dirac oscillator evolution operator 

We can understand the interferometric process clearer 
in the three-step process of fig. IH Here, the term e'^"'*^^ 
represents a beam splitter at the entrance of the interfer- 
ometer 



cos Vni 
-i sin dr. 



(30) 



the following term c describes the dephasing 

process in the two arms of the interferometer 











(31) 



and the remaining term stands for the final beam 

splitter which produces the interference between the de- 
phased beams that have traveled through different paths 
of the interferometer 



cos i sin 
i sin cos 0„i 



(32) 



Remarkably enough, this three-step process captures the 
essence of the relativistic dynamical properties in the 
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two-dimensional Dirac oscillator. The two incoming 
beams might be interpreted as the upper and lower com- 
ponents of the relativistic spinor. The first beam split- 
ter is responsible for the mixture of these components 
so that the two arms of the interferometer can be as- 
sociated to positive- and negative-energy solutions ([6]). 
In such manner, the time evolution inside the interfer- 
ometer can be understood as a phase shift between the 
positive- and negative-energy solutions since their phases 
evolve with opposite sign. Finally, the second beam split- 
ter is responsible for the interference of the interferom- 
eter beams, which consequently represents the interfer- 
ence of positive- and negative-energy solutions. This is 
exactly the essence of the Zitterbewegung in Relativis- 
tic Quantum Mechanics which can be surprisingly 
identified with a simple interferometric mechanism for 
the two-dimensional Dirac oscillator. 
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FIG. 4: "(Color online)" Schematic Mach-Zehnder interfer- 
ometer diagram of the Dirac oscillator evolution operator 

In the non-rclativistic regime described in section [TTTl 
the interferometric picture is significantly simplified. In 
this limit, the parameter rjm satisfies 7]^^ <C 1 for signi- 
ficative values of the initial number of orbital quanta ni 
. Therefore, we can approximate 



1 + -n^ e" 



1 
1 



(33) 



and thus the rotation angle « becomes vanishingly 
small. In this manner, the action of the two beam split- 
ters in fig. [3] is negligible, and the non-trivial remaining 
effect is the dephasing of the upper and lower components 
of the relativistic spinor in Eq. , which yields 



(f>ni (t) 



2ujni 



(34) 



As discussed in the preceding section, these phase shifts 
can manifest themselves in a Ramsey interferometer pro- 
viding a practical realization of Zitterbewegung. Beyond 
this example that requires the preparation of the atom in 
a superposition state, we can obtain a further example 
of Zitterbewegung dynamically induced by a first-order 
relativistic correction. 



The remarkable advantage of the interferometric inter- 
pretation developed in this section is the possibility to go 
beyond this non-relativistic regime, and consider higher- 
order relativistic corrections to the aforementioned dy- 
namics. In order to do so, we retain the relativistic cor- 
rections up to the following order 



1 



1 _ Z^i i^n. 

S 2 , 

^2 ^ 8-1 



which leads to the following evolution operator 



(35) 



■ sm ( 



Vni sm ( 



H + 



smc 



■ smc 



(36) 

Note that the rotation angle is no longer negligible, and 
the action of the beam splitters becomes noticeable up 
to the perturbative order considered so far. The inter- 
ference between the positive- and negative-energy solu- 
tions appears as a pure relativistic effect, leading to the 
spin-orbit Zitterbewegung if we consider an initial state 
■.^\ni-l)\xi) 

{L,)t = -A^nihsin' n^,t - h{ni - 1) + O {{^m)^) , 
{Jz)t= h{\-ni), 

(37) 

which clearly coincide with the expansion on the small 
parameter ^n; ^ 1 of the dynamical evolution described 
in Eqs. ([7]). The oscillations in the angular momentum 
observables are therefore a direct consequence of the in- 
terference between positive- and negative-energy solu- 
tions introduced by the Mach-Zehnder beam splitters. 
The visibility of the interference phenomenon in the spin 
degrees of freedom is 



V 



i^m + O ii^nif) , (38) 



which clearly fulfils V ^ 1 at this level of perturbation 
theory. As we discuss in appendix [B] the visibility of 
these Zitterbewegung oscillations increases considerably 
as relativistic effects become more pronounced, and sub- 
sequent perturbative orders are taken into account. 

Additionally, one must also consider the difference in 
the order of magnitude of the superposed frequencies in 
Eq. ([34|) . where 2ujni <C mc^/h, which makes it diffi- 
cult to observe the aforementioned instantaneous oscilla- 
tions. In such case, we can also perform a time average 
of Eqs. dSZl) 



= ~2^nih~hini-l) + 0{i^niy') , 

-~ 2(nih~l+0{{^niy), (39) 

{Jz)t^ h{\-ni), 



{S. 



which can be readily interpreted as a frequency shift in 
a non-relativistic left-handed harmonic oscillator lo 
uj + 2i{ see fig.©. 
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VI. CONCLUSIONS 

In this paper wc have considered the intriguing rela- 
tivistic Zitterbewegung in the two-dimensional Dirac os- 
cihator from an interferometric point of view. The exact 
mapping between the relativistic model and the Anti- 
Jaynes-Cummings interaction, suggests the use of quan- 
tum optical tools in the study of relativistic quantum 
phenomena. In this sense, the non-relativistic limit of the 
Dirac oscillator can be understood as a Ramsey interfer- 
ometric effect, and interesting Zitterbewegung-dyneLmics 
arise when the initial state is carefully prepared. Actu- 
ally, we have described how Zitterbewegung-free circular 
orbits become elliptic trajectories due to the interference 
of positive- and negative-energy modes. This insightful 
interferometric interpretation can be carried further on 
to subsequent relativistic corrections, in terms of a Mach- 
Zchnder interferometric process as shown in appendix [B] 
The effect of the Mach-Zehnder beam splitters is respon- 
sible of the spin-orbit Zitterbewegung^ and becomes more 
relevant as the relativistic parameter is increased. 

It is interesting to point out that all these exotic rel- 
ativistic effects may be observed in an ion-trap tabletop 
experiment following the proposal described in [l^ . In 
this experimental setting, the relativistic parameter can 
attain all possible values regarding current technology 
possibilities. This fact should allow the experimentalist 
to study this non-relativistic regime and the interfero- 
metric effects discussed in this work. 

Finally, we would also like to stress that an exciting di- 
alogue between Quantum Optics and Relativistic Quan- 
tum Mechanics scientists can be performed in the light 
of our results. These two communities can collaborate 
in order to offer a different perspective to archetypical 
phenomena of both disciplines. 

Acknowledgements We acknowledge financial support 
from the Spanish MEG project FIS2006-04885, the 
project CAM-UCM/910758 (AB and MAMD) and the 
UCM project PRl-A/07-15378 (AL). Additionally, we 
acknowledge support from a FPU MEG grant (AB), 




-h{ni - 1) 



FIG. 5: "(Color online)" Effective frequency displacement of 
a non-relativistic harmonic oscillator due to first order rela- 
tivistic corrections 



and the ESF Science Programme INSTANS 2005-2010 
(MAMD). 



APPENDIX A: STANDARD NON-RELATIVISTIC 
LIMIT 

In this appendix we derive the non-relativistic limit 
of the two-dimensional Dirac oscillator using the stan- 
dard techniques in Relativistic Quantum Mechanics [l6j . 
Let us then consider equation ([3]), where the relativis- 
tic spinor can be rewritten as j^*) :— [\ipi), |'!/'2)]*e~''^*/''. 
Then, equation ^ becomes a set of coupled equations 

(E - mc^)\ipi) = c [{px + imujx) - i{py + irriLuy)] \ip2), 
(E + TOC^)|V'2> = c [{px - imuJx) + i{py - imujy)] \^pl). 

(Al) 

In order to obtain the non-relativistic limit, we must 
derive the associated Klein-Gordon equations, which by 
virtue of the canonical commutation relations [a:;j,pfc] 
ihSjk, and using the definition of the orbital angular mo- 
mentum operator '■= xpy — ypx, become 

[E^ - m^c^)\i^{) = 2mc2 [H^^ - hoj - ojL,] 

[E^ - m^c'')\^2) = 2mc2 [iJ^D _^ _ ^^^j j^^)^ 

where we immediately identify the two-dimensional 
isotropic harmonic oscillator H^^ = /2m + mtj^r^/2. 
This fact already shows the connection between the 
relativistic Dirac oscillator in Eq. ([3]) and the usual 
non-relativistic harmonic oscillator. The non-relativistic 
regime is attained when the relevant energies in the sys- 
tem are negligible in comparison with the rest mass en- 
ergy. For the IV'i) component, we let E = mc^ + e where 
e ^ mc^, so that 

{E^ - m^c^) « 2mc^e + O (e^) . (A3) 

Substituting Eq. (|A3p in the Klein-Gordon equa- 
tion (|A2p . we obtain the corresponding non-relativistic 
limit of the two-dimensional Dirac oscillator 

e\^i)^[Hf^ -tiLj-uL^W^^). (A4) 

Finally, recovering the total energy E = mc^ -I- e, we 
obtain the following effective Hamiltonian for the non- 
relativistic limit 

H^^ - mc^ + [Hl^ -hio- uL,] , (A5) 

where the two-dimensional harmonic oscillator appears 
together with the orbital angular momentum. This pro- 
cedure must also be applied to the lower component 
I ■02), where the non-relativistic limit is attained setting 
E = —mc^ + e where e <IC mc^ 

{E^ - m^c^) w -2mc^t + O (e^) . (A6) 
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Substituting once more Eq. (|A6|) in the corresponding 
Klein-Gordon equation (|A2p . we obtain 



(A7) 



which directly leads to the effective Hamiltonian in the 

non-relativistic limit by restoring the original energy E ~ 

2 I 



H. 



(1) 



-mc 



(A8) 



where the usual two-dimensional harmonic oscillator 
arises naturally in this non-relativistic regime. Finally, 
using the chiral operators in Eq. ([4]), the two-dimensional 
harmonic oscillator can also be expressed in terms of 
these operators 



(alar + a\ai + , (A9) 



which leads to the corresponding non-relativistic effective 
Hamiltonian 



mc^ + 2hjja\ai 






-mc^ — 2fiuj{ajai + 1) 



(AlO) 

By virtue of the commutation relations [ai,a|] = 1. we 
can rewrite Eq. (jAlOp as 



mc^ + 2fiujalai 





—mc^ — 2fiwaia] 



(All) 



which coincides with the previous derivation using quan- 
tum optical tools ([TO]) . In this sense, the insightful Quan- 
tum Optics perspective introduced in [iJl offers a better 
understanding of the non-relativistic limit which is con- 
densed in fig.[ll In this regime, the spinorial levels can be 
only coupled through virtual transitions which is trans- 
lated into a displacement of the energies (fTT)) . 

We finally want to note that the inclusion of rest mass 
energy terms ±7tic^ in Eq. (jAlip is necessary if one wants 
to treat both spinorial components simultaneously. In 
this limit, these components are associated to positive- 
and negative-energy solutions, and their simultaneous 
treatment as considered above is a relativistic effect. 



APPENDIX B: PERTURB ATIVE SERIES OF 
THE NON-RELATIVISTIC LIMIT 

In this appendix we derive the complete perturbative 
series that arises naturally in the non-relativistic limit. 
Remarkably, we can obtain the aforementioned perturba- 
tive expansion to every order O ((C?^;)*') ^'^'^ '"^ Phys- 
ical interpretation of the corresponding terms. It turns 
out that the whole perturbative series can be interpreted 
in terms of dynamical Stark shift terms introduced in 
Sect. Illli and interferometric Ramsey processes as those 
discussed in Sect. |Vl In this sense, a complete descrip- 
tion of the phenomenology in the non-relativistic regime 



of the Dirac oscillator can be accomplished by only con- 
sidering the first two corrections in Sects. [1111 andlVl In 
light of these results, we claim that Zitterbewegung-like 
processes of the Dirac oscillator to any order can be fully 
described with the results of this work. 

Let us consider the unitary time evolution operator in 
Eq. ((28|) . which can be readily expressed as 

Uni = C/^"'' + 2i sin (t) sin 6'„, {sm9ni<Jz - cos 6l„,cry) , 

(Bl) 

where := e represents the zeroth order time 

evolution corresponding to the non-relativistic limit dis- 
cussed in Sect. mil This time evolution can be interpreted 
as the dephasing process inside the two arms of the Mach- 
Zehndcr interferometer ( see fig. HI). The remaining term 
contains the whole perturbative series and therefore we 
must expand in powers of the small parameter ^ni <^ 1 
to obtain the different corrections to the non-relativistic 
limit. Considering the following expansions 



sm' 



sin0„, cos0„, - i^niy/^ ~ 2{(nif/^ + 6(^nz)'/'- 
-20i^n,)y^+o[i^n,f/'), 

(B2) 

the time evolution operator (|Bip can be expressed as fol- 
lows 

C/„, = C/W + 2isin0„,(t) [-(^nO'/Vy + iini)a,+ 
+ 2{^nif/-'<7y - 3{^ni)^<7, - 6{^nif/^ay + 



O 



(B3) 



We observe from this expression how the subsequent per- 
turbative corrections present either a cTz term associated 
to a dispersive Stark-Shift dynamics, or a Uy term which 
can be identified with a Ramsey interference effect where 
the spinorial components get dynamically mixed dur- 
ing unitary evolution. Therefore, the perturbative series 
might be represented diagrammatically as in fig. |6l 

As we mentioned at the beginning of this appendix, the 
complete perturbative series can be interpreted in terms 
of the dynamical Stark shift term which accounts for the 
non-relativistic limit in section|TTTl and the Ramsey inter- 
ference term which accounts for the first order correction 
in section |V| Therefore, the full phenomenology and the 
characterization of Zitterbewegung can be accomplished 
regarding these two regimes. 

Finally, let us interpret the different terms of the so- 
called perturbative expansion with the language of Quan- 
tum Optics. The unitary operator in Eq. (|B3[) describes 
the time evolution inside the invariant subspace 7i„, in 
Eq. (|13|1 . If we recover the full Hilbert space descrip- 
tion, we obtain that the even terms of the perturbative 
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expansion are 



2k 



\2k sin (pm (t) 











(B4) 



where k = 0,1... and we have mtroduced certain cou- 
phng constants A2fc(t) := C2k£,^ which involve increasing 
powers of the relativistic parameter, and C2k G C follow 
directly from Eqs. ()Bip . ()B3p . These terms represent a 
kind of dynamical Stark shift between the spinorial lev- 
els proportional to the fc— th power of the orbital quanta 
number n/. They might be considered as certain shifts 
produced by 2A;-virtual transitions between the spinorial 
levels. 

In the same manner, we can rewrite the odd terms as 



u- 



2k+l 



X2k+i sin 4>niit) 



(a|aO'a| 
-aiiajai)'' 



where the coupling parameter X2k+i — C2k+i^^~ be- 
comes more important as the relarivistic parameter in- 
creases, and C2fc+i G C also follow from Eqs. (|B1[) . (|B3|) . 
These pcrturbative terms (jB5l) can be directly expresses 
as a generalized Anti-Jaynes-Cummings evolution 



U 



2k+l 



A 



2k+l 



-A, 



Ansin0„,(t), (B6) 



where we have introduced the bosonic operator 

Ai, := aiiajai)''. (B7) 



This effective time evolution (|B6p can be interpreted as a 
generalized AJC-model where intensity of the couplings 
depends on the fc— th power of the number of left-handed 
quanta n/ [H, 




FIG. 6; "(Color online)" Diagrammatic scheme of the non- 
relativistic perturbative series. The even terms represented by 
two parallel lines correspond to the az coupling associated to 
a dynamical Stark shift term, whereas the odd terms depicted 
by crossing lines stand for the ay interaction associated to a 
Ramsey interferometric effect. Note that the different scales 
stress the relevance of each term in the non-relativistic regime. 
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